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I review the phenomenological and cosmological constraints on the parameter 
space associated with the Constrained Minimal Supersymmetric Standard Model 
(CMSSM). The effect of the recent WMAP determination of the cold dark matter 
density will be discussed. A very constrained model (based on minimal supergrav- 
ity relations between the bi- and tri-linear supersymmetry breaking mass terms) 
will be outlined. 



1 Introduction 

It is well known that supersymmetric models with conserved i?-parity contain 
one new stable particle which is a candidate for cold dark matter 
There are very strong constraints, however, forbidding the existence of stable 
or long lived particles which are not color and electrically neutral. Strong 
and electromagnetically interacting LSPs would become bound with normal 
matter forming anomalously heavy isotopes. Indeed, there are very stronj 
upper limits on the abundances, relative to hydrogen, of nuclear isotope; 
n/n H < 10- 15 to 1(T 29 for 1 GeV < m < 1 TeV. There are relatively few 
supersymmet ric c andidates which are not colored and are electrically neutral. 
The sneutrind^G] j s orLe possibility, but in the MSSM, it has been excluded 
as a dark matter candidate by direct!^ and indirect searches. Another pos- 
sibility is the gravitino which is probably the most difficult to exclude. I will 
concentrate on the remaining possibility in the MSSM, namely the neutrali- 
nos. 

There are four neutralinos, each of which is a linear combination of the 
R = — 1, neutral fermionffl the wino W 3 , the partner of the 3rd component 
of the SU(2)l gauge boson; the bino, B, the partner of the U(l)y gauge 
boson; and the two neutral Higgsinos, Hi and H^- In general, neutralinos can 
be expressed as a linear combination 

X = aB + (3W 3 + 7#i + SH 2 (1) 

The solution for the coefficients a, /3, 7 and S for neutralinos that make up the 
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LSP can be found by diagonalizing the mass matrix 
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where Mi(M2) are the soft supersymmetry breaking U(l) (SU(2)) gaugino 
mass terms, fx is the supersymmetric Higgs mixing mass parameter and since 
there are two Higgs doublets in the MSSM, there are two vacuum expectation 
values, vi and v 2 - One combination of these is related to the Z mass, and 
therefore is not a free parameter, while the other combination, the ratio of 
the two vevs, tan/3, is free. 

The most general version of the MSSM, despite its minimality in par- 
ticles and interactions contains well over a hundred new parameters. The 
study of such a model would be untenable were it not for some (well moti- 
vated) assumptions. These have to do with the parameters associated with 
supersymmetry breaking. It is often assumed that, at some unification scale, 
all of the gaugino masses receive a common mass, m 1 / 2 - The gaugino masses 
at the weak scale are determined by running a set of renormalization group 
equations. Similarly, one often assumes that all scalars receive a common 
mass, mo, at the GUT scale (though one may wish to make an exception 
for the Higgs soft masses). These too are run down to the weak scale. The 
remaining supersymmetry breaking parameters are the trilinear mass terms, 
Ao, which I will also assume are unified at the GUT scale, and the bilinear 
mass term B. There arc, in addition, two physical CP violating phases which 
will not be considered here. 

The natural boundary conditions at the GUT scale for the MSSM would 
include fx, the two soft Higgs masses (mi and mi) and B in addition to mi/ 2 , 
mo, and Aq. In this case, upon running the RGEs down to a low energy 
scale and minimizing the Higgs potential, one would predict the values of Mz, 
tan/3, and the Higgs pseudoscalar mass, tua (in addition to all of the sparticle 
masses). Since Mz is known, it is more useful to analyze supersymmetric 
models where Mz is input rather than output. It is also common to treat 
tan/3 as an input parameter. This can be done at the expense of shifting 
fi (up to a sign) and B from inputs to outputs. When the supersymmetry 
breaking Higgs soft masses are also unified at the GUT scale (and take the 
common value mo), the model is often referred to as the constrained MSSM 
or CMSSM. Once these parameters are set, the entire spectrum of sparticle 
masses at the weak scale can be calculated. 
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2 Constraints 



2.1 The Relic Density 

The relic abundance of LSP's is determined by solving the Boltzmann equation 
for the LSP number density in an expanding Universe. The technique^ used is 
similar to that for computing the relic abundance of massive neutrinos^. The 
relic density depends on additional parameters in the MSSM beyond m^, /A 
and tan/3. These include the sfermion masses, wif, as well as m^, all derived 
from mo, A®, and mm. To determine the relic density it is necessary to obtain 
the general annihilation cross-section for neutralinos. In much of the param- 
eter space of interest, the LSP is a bino and the annihilation proceeds mainly 
through sfermion exchange. Because of the p-wave suppression associated 
with Majorana fermions, the s-wave part of the annihilation cross-section is 
suppressed by the outgoing fermion masses. This means that it is necessary to 
expand the cross-section to include p-wave corrections which can be expressed 
as a term proportional to the temperature if neutralinos are in equilibrium. 
The partial wave expansion results in a very good approximation to the relic 
density expect near s-channel annihilation poles, thresholds and in regions 
where the LSP is nearly degenerate with the next lightest supersymmetric 
particle^, where a more accurate treatment is necessary. 

The preferred range of the relic LSP density has always been restricted 
to a relatively narrow range 0.1 < JIcdm^ 2 < 0.3, where values much smaller 
than the lower bound are disfavoured by analyses of structure formation in the 
CDM framework and the upper limit assumes only that the age of the Universe 
is > 12 Gyr. However, one should note that the LSP may not constitute all 
the CDM, in which case Olsp could be reduced below this value. This range 
has, however, been altered significantly by the recent improved determination 
of the allowable range of the cold dark matter density obtained by combining 
WMAP and other cosmological data: 0.094 < Qcdm < 0.129 at the 2-cr 
leveP. 



2.2 Accelerator Bounds 

The most relevan t co nstraints on the su persymmetric parameter space are: 
m x ± > 104 GevEl, m s > 99 Gevldand m h > 114 GeV EH. The for- 
mer two constrain 111^2 and mo directly via the sparticle masses, and the 
latter indirectly via the sensitivity of radiative corrections to the Higgs m ass 
to the sparticle masses, principally mrr. The latest version of FeynHiggsEl 
is used for the calculation of m^. The Higgs limit imposes important con- 
straints principally on particularly at low tan/3. Another constraint is 
the requirement that the bran ching ratio for b — ► sj is consistent with the 
experimental measurements^! These measurements agree with the Standard 
Model, and therefore provide bounds on MSSM particles, such as the chargino 



3 



and charged Higgs masses, in particular. Typically, the b — > S7 constraint is 
more important for /1 < 0, but it is also relevant for fi > 0, particularly 
when tan/3 is large. The constraint imposed by measurements of b — > sj also 
excludes small values of mm- 

Finally, there are regions of the (mi/2, mo) plane that are favoured by the 

BNL measurement^! of g^ — 2 at the 2-a level, corresponding to a deviation of 
(33. 9± 11.2) x 1CT 10 from the Standard Model calculatiorP^I using e + e ~ data. 
One should be however aware that this constraint is still under discussion and 
it will not be used to constrain tan (3. All the fi > planes would be consistent 
with — 2 at the 3-u level, whereas /x < is disfavoured even if one takes a 
relaxed view of the g^ — 2 constraint. 



3 Results for the TO1/2 — mo Plane Before and After WMAP 

As noter earlier, the spectrum for a given model in the CMSSM is determined 
by the parameter set mi/2, mo, tan/3, Aq, and the sign of fj,. For now, I will 
assume that A = 0. Then, for a given value of tan/3 and sgn(/i), the resulting 
constraints can be displayed on the TO]y 2 — too plane. Let us first consider the 
region of parameter space with relatively low values of m 1 / 2 and mo, i.e. the 
'bulk' region in the CMSSM for tan/3 = 10 and /j, > shown in Fig. [J 18 . 
The light shaded region correspond to 0.1 < Q, x h 2 < 0.3. The dark shaded 
region has to— < m x and is excluded. The light dashed contours indicate 
the corresponding region in Q x h 2 if one ignores the effect of coannihilations. 
Neglecting coannihilations, one would find an upper bound of ~ 450 GeV on 
my 2 , corresponding to an upper bound of roughly 200 GeV on trig. 

Coannihilations are important when there are several particle species i, 
with different masses, and each with its own number density rii and equilib- 
rium number density no,i- hi this case inclusion of coannihilations is required 
because of the degeneracy of \ and t 1 . The effect of coannihilations is to 
create an allowed band about 25-50 GeV wide in to for my 2 < 1400 GeV, 
which tracks above the m^ = m x contoujCSCl. 

A larg er view of th e tan/3 = 10 parameter plane is shown in the left panel 
of Fig. [2|| 18 | 20 | 21 | 22 | i ncm d ec [ nere are the phenomenological constraints 
discussed above. As one can see, the region preferred by g — 2 overlaps very 
nicely with the 'bulk' region for tan /3 = 10 and fi > 0. In the second panel 
of Fig. |21 we see the effect of imposing the WMAP range on the neutralino 
densitji^lMI. We see immediately that (i) the cosmological regions are gen- 
erally much narrower, and (ii) the 'bulk' regions at small rowj and too have 
almost disappeared, in particular when the laboratory constraints are im- 
posed. Looking more closely at the coannihilation regions, we see that (iii) 
they are significantly truncated as well as becoming much narrower, since the 
reduced upper bound on Q, x h 2 moves the tip where m x = n%f to smaller TO4/2 
so that the upper limit is now m.1/2 Ss 950 GeV or m x < 400 GeV. 
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Figure 1. The light-shaded 'bulk' area is the cosmologically preferred region with 
0.1 < Q x h 2 < 0.3. The light dashed lines show the location of the cosmologically preferred 
region if one ignores coannihilations with the light sleptons. In the dark shaded region in 
the bottom right, the LSP is the f\ , leading to an unacceptable abundance of charged dark 
matter. Also shown is the isomass contour m x ± = 104 GeV and = 110, 114 GeV, as 
well as an indication of the slepton bound from LEP. 



Another mechanism for extending the allowed CMSSM regio n to lar ge m x 
is rapid annihilation via a direct-channel pole when m x ~ \rriA 1 ^ 1 ^ H Since 
the heavy scalar and pseudoscalar Higgs masses decrease as tan (3 increases, 
eventually 2m x ~ jtia yielding a 'funnel' extending to large m 1 ^ 2 and too at 
large tan (3, as seen in Fig. The difficulty and necessary care involved in 
calculations at larg e tan/3 were previously discussed^. For related CMSSM 
calculations see 1 ^". 

In the second panel of Fig. EJ we see the effect of imposing the WMAP 
range on the neutralino density^Sl We see rapid-annihilation funnels that are 
also narrower and extend to lower toi / 2 and Too than previously. They weaken 
significantly the upper bound on m x for tan/3 > 50 for fi > 0. 

Shown in Fig. 0] are the wmap 

lineal 

of the (m 1 /2,mo) plane allowed 
by the new cosmological constraint 0.094 < fl x h 2 < 0.129 and the labora- 
tory constraints listed above, for fi > and values of tan (3 from 5 to 55, in 
steps A(tan/3) — 5. We notice immediately that the strips are considerably 
narrower than the spacing between them, though any intermediate point in 
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Figure 2. The (wii/2i m o) planes for (a) tan j3 = 10 and /i > 0, assuming Ao = 0,mt = 
175 GeV and mt(mi)^ = 4.25 GeV. The near-vertical (red) dot-dashed lines are the 
contours = 114 GeV as calculated using FeynHiggs and the near-vertical (black) 
dashed line is the contour rn x ± = 104 GeV. Also shown by the dot-dashed curve in the 
lower left is the corner excluded by the LEP bound of mg > 99 GeV. The medium (dark 
green) shaded region is excluded by b — > s'y, and the light (turquoise) shaded area is the 
cosmologically preferred regions with 0.1 < fl x h 2 < 0.3. In the dark (brick red) shaded 
region, the LSP is the charged f\ . The region allowed by the E821 measurement of a M at 
the 2-a level, is shaded (pink) and bounded by solid black lines, with dashed lines indicating 
the 1-a ranges. In (b), the relic density is restricted to the range 0.094 < Q. x h 2 < 0.129. 
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Figure 3. As in Fig. \^for tan /3 = 50 and fj, > 
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Figure 4. Tfte strips display the regions of the ("H/2, "io) plane that are compatible 

1,2 



5, 10, 15, 20, 25, 30, 35, 40, 45, 50, 55. 
2-o~ level have darker shading. 



The parts of the strips compatible with cj M — 2 at the 



the (?Bi/2, ?rio) plane would be compatible with some intermediate value of 
tan/3. The right (left) ends of the strips correspond to the maximal (minimal) 
allowed values of m 1 / 2 and hence m x . The lower bounds on m 1( / 2 are due to 
the Higgs mass constraint for tan/3 < 23, but are determined by the b — ^ sj 
constraint for higher values of tan/3. The upper bound on miy 2 for tan/3 > 50 
is clearly weaker, because of the rapid-annihilation regions. 



4 Using SUGRA to constrain the CMSSM 

While the CMSSM models described above are certainly mSUGRA inspired, 
minimal supergravity models can be argued to be still more predictive or 
constrained. Let us begin by assuming that supersymmetry is broken in a 
hidden sector so that the superpotential can be written as a sum of two terms, 
W ~ F(<p) + g(C), where <j> represents all observable fields and ( all hidden 
sector fields. Then the scalar potential in an N = 1 supergravity theory with 
minimal Kahler potential K((j>, () = <fi<f)* + ((* + In | W\ 2 , can be written a 
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+ <f>*(9(0 + F{m 2 - 3|g(C) + F(cj>)f 



(3) 



We furthermore must choose g(() such that when £ picks up a vacuum ex- 
pectation value, supersymmetry is broken. When the potential is exp anded 
and terms inversely proportional to Planck mass are dropped, one findJ^EHI 
1) scalar mass universality with too = (g); 2) trilinear mass universality with 
A = (dg/dQ(0 + (g)(C) 2 ; and 3) B = A - m . 

In fact, one of the primary motivations for the CMSSM, and for scalar 
mass universality in p art icular, comes from the simplest model for local su- 
persymmetry breakinj^ill which involves just one additional chiral multiplet 
£ with a superpotential of the form 

g(0 = KC + /3) (4) 

with |/3| = 2 - a/3, ensuring that the cosmological constant A = 0. After 
insertin g th e vev for £, (£) = Vo — 1. The scalar potential in this model takes 
the fornffl 



y _ e (4-2V5) 
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\v+{VZ-l){v + F{^))\ 
+ \^- + ^{v + F{4>))\ 2 -Z\v + Fm 2 



e (4-2v / 3)|^|2 



+TO 3/2 e( 2 -^ (0— - V3F + h.c.)) + m 2 3/2 W* , (5) 



where the gravitino mass is given by to 3 / 2 = e 2 ^v. 

First, up to an overall rescaling of the superpotential, F — * e^~ 2 F, the 
first term is the ordinary i^-term part of the scalar potential of global super- 
symmetry. The next term, which is proportional to m 3 / 2 , provides universal 
trilinear soft supersymmetry- breaking terms A = (3 — v / 3)to 3 / 2 and bilinear 
soft supersymmetry- breaking terms B — (2 — V / 3)'Ti3/2j i-e., a special case of 
the general relation above between B and A. Finally, the last term represents 
a universal scalar mass of the type advocated in the CMSSM, with to 2 , = w 2 ^ 2 , 
since the cosmological constant A vanishes in this model, by construction. 

Given a relation between Bq and Aq, we can no longer use the standard 
CMSSM boundary conditions, in which mi/ 2 , too, Aq, tan/3, and sgn(n) 
are input at the GUT scale with fi and B determined by the electroweak 
symmetry breaking condition. Now, one is forced to input Bn and instead 
tan/3 is calculated from the minimization of the Higgs potential" . It should 
be stressed that this type of model is one that truly originates from minimal 
supergravity. 

In Fig. [31 the contours of tan/3 (solid blue lines) in the (toi/ 2 , too) planes 
for two values of A = Aq/toq, B — Bo/mo = A — 1 and the sign of /i are 
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Figure 5. Examples of (mj/j , mo) planes with contours of tan /3 superposed, for /i > and 
(a) the simplest Polonyi model with A = 3 — \/3, B = A — 1 and (b) A = 2.0, B = A — 1. 
7n eac/i panel, we show the regio ns e xcluded by the LEP lower limits on MSSM particles, 
those ruled out by b —> s"/ decay^^ (medium green shading), and those excluded because 
the LSP would be charged (dark red shading). The region favoured by the WMAP range 
flCDAih 2 = 0-H26^Q'QQgj has light turquoise shading. The region suggested by g^ — 2 is 
medium (pink) shaded. 



displayed*"*!. Also shown are the contours where m x ± > 104 GeV (near- 
vertical black dashed lines) and rrih > 114 GeV (diagonal red dash-dotted 
lines). The excluded regions where m x > have dark (red) shading, those 
excluded by b — ► 57 have medium (green) shading, and those where the relic 
density of neutralinos lies within the WMAP range 0.094 < tl x h 2 < 0.129 
have light (turquoise) shading. Finally, the regions favoured by — 2 at the 
2-cr level are medium (pink) shaded. 

In panel (a) of Fig. we see that the Higgs constraint combined with 
the relic density requires tan/3 > 11, whilst the relic density also enforces 
tan/3 < 20. For a given point in the mi/2 — wo plane, the calculated value 
of tan /3 increases as A increases. This is seen in panel (b) of Fig. when 
A = 2.0, close to its maximal value for fi > 0, the tan/3 contours turn over 
towards smaller m.1/2, and only relatively large values 25 < tan/3 < 35 are 
allowed by the b — > 57 and flcDMh 2 constraints, respectively. For any given 
value of A, there is only a relatively narrow range allowed for tan /3. 

At higher values of A (and high tan/3), the off-diagonal elements in the 
squark mass matrix become large at large mo. Therefore, we find no solutions 
which are phenomenologically viable above A ~ 2.5. This is because the 
regions where the LSP is the r or the t close off the parameter space. In fact, 
this feature is generic in the CMSSM as shown in Fig. 3 of ^3 This effect is 
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more severe at large tan (3, which further compounds the difficulty in going to 
large values of A in the type of models discussed here. At very low values of 
A (< —1.9), there are no viable solutions due to the competition between the 
Higgs mass bound and the relic density. At low A, solutions give low values 
of tan/3 which in turn give low Higgs masses unless m 1 / 2 is very large. For 
A < —1.9, m 1 / 2 is pushed passed the endpoint of the coannihilation region. 

In addition, we note the absence of the funnel regions. This is due to the 
relatively small values of tan j3 allowed in the class of models considered here: 
we recall that the funnel region appears only for large tan /3 > 45 for /i > 
and tan (3 > 30 for [i < in the CMSSM. 

Rather than further constrain the CMSSM as described above , one 
can generalize the CMSSM case to include non-universal Higgs masse d I I 
(NUHM), in which case the input parameters include /j, and m^, in addition 
to the standard CMSSM inputs. In this case, the two soft Higgs masses, 
mi, W2 are no longer set equal to tjiq and are calculated from the electroweak 
symmetry breaking conditions. The NUHM paramete r sp ace was recently 
analyzed^ and a sample of the results were reviewed in l^l 
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